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Abstract 

A Hamiltonian cycle of a graph is a closed path that visits every vertex 
once and only once. It has been difficult to count the number of Hamiltonian 
cycles on regular lattices with periodic boundary conditions, e.g. lattices on a 
torus, due to the presence of winding modes. In this paper, the exact number 
of Hamiltonian cycles on a random trivalent fat graph drawn faithfully on a 
torus is obtained. This result is further extended to the case of random graphs 
drawn on surfaces of an arbitrary genus. The conformational exponent 7 is 
found to depend on the genus linearly. 
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1 Introduction 



Properties of Hamiltonian cycles and walks on regular and random lattices 
have been attracting much attention recently [1, 2, 3, 4, 5, 6]. A Hamiltonian 
cycle (walk) of a graph is a closed (open) path which visits every vertex once 
and only once, i.e., a self-avoiding loop (walk) which visits all the vertices. 
The system of a Hamiltonian cycle (walk) on a lattice serves as a model of 
compact ring (linear) polymer which fills the lattice completely [7]. It is also 
relevant to the protein folding problem [8, 9, 10]. 

Among the properties of Hamiltonian cycles, the number of Hamiltonian 
cycles on a given graph is one of the most fundamental and interesting quan- 
tities. The number is directly related to the entropy of a lattice polymer in 
the compact phase. When it is non-zero, it also coincides with the degen- 
eracy of optimal solutions to the traveling salesman problem on the graph. 
The number has been calculated exactly or numerically for a number of fixed 
regular lattices [11, 12, 13, 14, 15, 16, 1]. 

In this article, I am concerned with random lattices and stTidy the number 
of Hamiltonian cycles on them. Specifically, I exactly evaluate the number 
^Q^s^,g 'H{G) / {^Aut G), where the number of Hamiltonian cycles 
on a given graph G is denoted by H{G). The ensemble S'^'^ is the set of all 
trivalent fat graphs that have n vertices and can be drawn faithfully on a 
surface of genus g ( but not on that of g — 1). Sec Fig. 1 for an example. 
The integer ^^Aut G is related to the symmetry of G and is defined precisely 
in section 2. This result extends those of refs. [5, 6] where only the planar 
case {g — 0) has been studied. 

The case g — 1 (torus) is especially interesting in that it can be compared 
with the corresponding problem on a two-dimensional fixed regular lattice 
with the periodic boundary condition^ for each of the two directions [2]. This 
comparison provides a good illustration of nature of statistical systems on 
random lattices. 

Prom the exact integral expression of Fn^\ one can extract the large- 
n asymptotics of the 'random average' of 1-L{G). It is found that the site 
entropy is independent of the genus g while the conformational exponent 
7 depends on it linearly. This behavior is consistent with the KPZ-DDK 
scaling [17, 18, 19] for two-dimensional gravity coupled to c = — 2 conformal 

^In some literatures, the loops is said to be a walk satisfying the periodic boundary 
condition. In the present context, however, it is imposed on lattices in the absence of loops 
or walks. 
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Figure 1: An example of a trivalent fat graph G G 5''^'^ which can be drawn 
on a torus but not on a sphere. A Hamiltonian cycle is drawn in a gray line. 
For this graph, H{G) — 8 because a direction (arrow) and a base point (dot) 
are associated with a Hamiltonian cycle. 



matter [5, 6]. 

The organization of the paper is as follows. I give definitions and fix 
notations in section 2. In section 3, the calculation for a random graph on a 
torus is presented, putting emphasis on the difference from the case of a fixed 
regular lattice on a torus. The analysis is extended to the case of surfaces of 
an arbitrary genus to yield a simple integral expression in section 4. I discuss 
my results in section 5. 

2 Definitions 

Definitions in this section generalize those given in ref. [6]. 

Let S"''^ be the set of all trivalent fat graphs that have n vertices pos- 
sibly with multiple-edges and self-loops and can be drawn faithfully on an 
orientable surface of genus g > (but not on S^-i). An example of 
a trivalent fat graph is drawn in Fig 1. Graphs that are isomorphic are 
identified. The set S'^'^ is the labeled version of 5""'^, namely, vertices of 
G e S'^'^ are labeled 1, . . . ,n and Gi, 6*2 G S"-'^ are considered identical only 
if a graph isomorphism preserves labels. The symmetric group of degree n 
naturally acts on S"''^ by the label permutation. The stabilizer subgroup of 
G is called the automorphism group Aut G. 

A Hamiltonian cycle of a labeled graph G e S'^'^ is a directed closed path 
(consecutive distinct edges connected at vertices) which visits every vertex 
once and only once. Hamiltonian cycles are understood as furnished with 
a direction and a base point (denoted by an arrow and a dot in figures). 



2 



The number of Hamiltonian cycles of G is denoted by Ti.{G) because it is 
independent of the way of labehng. 
The quantity I study in this work is 



and the generating function 

oo 

F^^){p):^Y,p-F^^\ (2) 



n=0 



The key observation made in ref. [6] has been that Fn'' can be written as the 
number of isomorphism classes of the pair (graph, Hamiltonian cycle) and 
thus is an integer. It is also true for Fn"^ with g' > : 

F^s) ^ #({(6-, C"^)!^ e ^"'^ : Hamiltonian cycle on G}/ (3) 

where (Gi, C{) ~ ((52, C2 ) if and only if Gi and G2 are isomorphic (forgetting 
the labels) and the isomorphism maps onto with the direction and 
the base point preserved. 



3 Lattices on a torus 

Before proceeding to the calculation of Fn \ let us take a glance at the 
case of fixed regular lattices on a torus Si. Given a fixed base point, a 
Hamiltonian cycle represents an element of the fundamental group 7r(Ei) ~ 
Za+Z6 where a and h are the meridian and the longitudinal cycles. Reflecting 
this fact, there are many topological sectors for self- avoiding loops (thus 
for Hamiltonian cycles) as depicted in Fig 2. More precisely, an element 
mitt + 777.26 e 7ri(Si) can be represented by a self-avoiding loop if and only if 
the pair (7771,7712) satisfies 

(mi,7772) = (0,0),(±1,0),(0,±1), or gcdK, 7772] = 1. (4) 

This rich structure gives rise to an interesting question: how many Hamilto- 
nian cycles belong to a given topological sector ? But at the same time, the 
existence of many topological modes is an obstacle to calculating 7i{G). 
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Figure 2: Examples of the ways how self-avoiding loops (gray lines) wind 
around tori. A torus Ei is drawn as a rectangle with two pairs of edges 
identified, with the generator a, 6 G 7ri(Ei). From left to right, the loop 
represents 0, a, b, 2a + b, and 2a + 36 e 7ri(Ei), respectively 



In the transfer matrix method for calculating TC{G), one maps the system 
into a state sum model with a local weight. Once mapped to a state sum 
model, its transfer matrix can be diagonalized by employing the Bethe ansatz 
[14, 15] or numerical calculation [20, 16]. In the mapping, it is crucial to 
avoid contributions from vacuum loops, or small loops disconnected from 
the largest loop component in question. This non-local constraint can be 
imposed as follows. First one assigns directions to each loop component and 
makes the weight direction dependent. The sum over all the assignments is 
taken. Then the weight is tuned so as that the cancelation occurs between 
two configurations having the vacuum loop of the same shape but with the 
opposite directions. 

This procedure is straightforward for lattices with disk topology because 
all vacuum loops are topologically trivial. In the case of cylinder geometry, 
vacuum loops winding around the cylinder should be taken care of. To make 
the cancelation work, one introduces the 'seam' and associates additional 
weights to vacuum loops going across the seam and winding around the 
cylinder. 

When one further imposes periodicity to the other direction to have a 
torus, one come to have a huge number of topologically inequivalent vacuum 
loops as shown in Fig. 2. As far as the present author knows, there is 
no simple way to avoid a contribution from each of these vacuum loops in 
the transfer matrix method. Thus, for fixed lattices on a torus, one has to 
employ the direct enumeration or the field theoretic approximation in order 
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to evaluate 7i{G) [2]. 

Note that the Coulomb gas method [4, 1, 3], which is capable of deter- 
mining exact critical exponents, again works for the disk or cylinder topology 
but not for the torus topology ^. 

For random lattices, however, the situation changes dramatically. Chang- 
ing the order of the summation, I shall sum over random lattice structures 
first fixing the topological sector of the Hamiltonian cycle. Then I sum over 
the topological sectors. Because the summation is performed over random 
lattices, there is no preferred basis (a, 6). One is free to perform a modular 
transformation on a torus to map mia -\- 777,26 1— > 1 • a, 6 1— > 6 if mi 7^ or 
m2 7^ 0. In other words, only two topological classes of cycles are distin- 
guishable: the non-winding sector (mi = m2 = 0) and the winding sector 
(mi ^ or m2 7^ 0). A contribution from each sector is denoted by Fn''''^ 
and Fn^'^ ^ respectively (Fig 3) 

00 

FW = F«'0 + F«'\ fW'-(p):=5^p-F«'-. (5) 

n=0 

Imagine that one walks along a Hamiltonian cycle in the specified direc- 
tion starting from the base point and records the order of right and left turns. 
Then one can associate a sequence of T-shaped objects with the order as de- 
picted in Fig. 3. If there are k right turns and n — k left ones (0 < /c < n), 
there are (^) possible orderings. The Hamiltonian cycle goes through the 
horizontal segment of T while the vertical segment is left unvisited. The 
n vertical segments should be connected pairwisely to reproduce the graph 
completely. The pairs are connected with fat black lines on the given surface 
(Fig.3). 

Now I consider and F^'^'^ separately. 

For non- winding sector there should be no connection between the 

right hand side and the left hand side of the cycle because it divides the torus 
into a disk and that with a handle. Let A^^'* denote the number of ways of 
connecting k objects on a disk with g handles (but not on that with g — 1 

^In contrast to the case of the Hamiltonian cycles or the compact polymer, various 
critical exponents of the dense polymer on torus was exactly determined by Duplantier 
and Salcur [21]. This is because the dense polymer is universal in the sense that the 
exponents does not depend on details of the lattice. 
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Figure 3: Configurations contributing to each of two topoiogicai sectors (ieft: 
Fg^'*'^, riglit: Fq^^'^). The sequence of T's in tfie center is obtained by walking 
along the Hamiltonian cycle (the gray line in the T's). The vertical segments 
of T's on both sides should be connected by black fat hues on the surface. 



handles). I have 

On the other hand, for winding sector Fn^'^ I obtain 

fc=o ^ ^ 

where A^i^l^_f. is the number of ways of contracting k and n — k objects on each 
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end of a cylinder with g handles attached (and with at least one connection 

between two ends of the cylinder). 

It has been known that A's above can be written as connected correlation 
functions of the hermitian gaussian matrix model [22, 23]: 



V OO 

-trM'^-trMM =y N-^^-^A^']. 

/ c n—n 



(8) 
(9) 



c 9=0 

The expectation value (• • • ) is defined by 



(...):= I y dM(•••)e->^^ (10) 

where M is a,n N x N hermitian matrix variable and the normalization Z is 
fixed so as to satisfy (1) = 1. The subscript c in (• • • )c means the connected 
correlation function 

(tr/(M) tigiM)), := (tr/(M)tr<7(M)) - (tr/(M)) (tr(?(M)). (11) 

Because the measure in (10) is gaussian, Ai^^ and a[^1 can be obtained 
easily by the method of orthogonal polynomial [24, 25] or by the loop equation 
[26, 27]. It can be shown that 

oo 

^(0)(x) :=^x-'=-Mf = i (^p- V^^) , (12) 

fc=0 

oo 

A('\x) -J^x-'^-UW = (x^ -4)-5/^ (13) 



k=0 
oo 

k,e=o 



1 ^1-^) (l-^) A^ ^-^ l . (14) 



24 V V^^zriy I ^/y^^ J \ ^-y 
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Combining these expressions, I finally arrive at exact integral expressions for 



where the contours for x, y go around the cut [—2, 2] counterclockwise. Care- 
ful inspection of the singularities shows that 

F^^^'°(p) — ^^ — (| — p)~^ + (less singular terms), (17) 
F^^^'^{p) =—;^ — (t — p)^^ + (less singular terms). (18) 

2' • TT 

This implies that Pji^^''^ grows as 4"n° for n — > oo and that the winding and 
the non- winding sectors contribute with the ratio of 3 : 2. 

The combinatorial argument leading to (6) and (7) was first used by 
Duplantier and Kostov in analyzing the dense and dilute phases of polymers 
on a planar random lattice [28, 29]. Then it has been explicitly recognized 
that it can be equally applied to the compact polymer, or the Hamiltonian 
cycle problem on a planar random lattice [5, 6]. This the first time that the 
combinatorial argument is successfully applied to higher-genus cases. 



4 Lattices on surfaces of arbitrary genus 

For random lattices drawn on a surface Eg of genus g, I find 

^i" = ltt ("^ + 1 ± il) (19) 

h=0 k=0 ^ ^ fc=0 ^ ^ 

by enumerating all the topological configurations. Equation (19) can be 
understood as follows (See Fig. 4). When one cuts the surface along a 
Hamiltonian cycle, one is left with a (possibly disconnected) surface with two 
circular boundaries. If the surface remains connected, the surface obtained 
should be with two holes (the second term in Fig 4). If it becomes 

disconnected, one has E/^ and ^g-h which are to be glued together along the 
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Figure 4: Topologically inequivalent ways to wind a Hamiltonian cycle 
around a surface of genus g. 



Hamiltonian cycle to form (the first term in Fig 4). This exhausts all the 
possibihties. There are precisely [| + 2] topologically inequivalent sectors. 
Note that the summands for h = ho and h — g — ho in the first term in 
(19) belong to an identical topological sector. Equation (19) can be used to 
calculate contributions from each topological sectors as well as the sum Fri^ 
using (8) and (9). 

At this point, I am naturally led to define the 'all-genus' generating func- 
tion 

oo oo 

F(p,t) = ^i^^p"F(^). (20) 

g=0 n=0 

Notably, this simplifies the expression a lot. I find that 

^) = ^ Ep"E it) (^trM^^trM--) , (21) 

n=0 jk=0 ^ / \ / 

where N is the dimension of the matrix variable M. The correlation function 
in the right hand side includes both connected and disconnected parts. In 
fact, I can easily recover (19) by expanding (21) in 1/N. 
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I further rewrite (21) to have a simple integral expression. I obtain 




In the above derivation, I have made use of the integral expression of the 
F-function. It can be easily checked that the coefficients of A^~^ and in 
(22) exactly agree with eq. (3) and the result in ref. [6], respectively. 

The expression (22) is not only simple but also practically useful. It 
enables one to calculate Fn^^ by the method of orthogonal polynomials [24, 
25]. 

Because F{p,^) has been written as a correlation function of the ga.us- 
sian hermitian matrix model in (22), the system should fall into the same uni- 
versality class as the topological gravity or two-dimensional quantum gravity 
coupled to c = — 2 matter [30, Subsection 9.3]. Namely, the double scaling 
behavior 

np,w)-ii-pffiN'{\-pr) (23) 

or equivalently the large-n asymptotics 

F^«) ~ 4"n2^-2 {n ^ oo) (24) 

is implied. 

The central charge c = —2 was previously obtained by Duplantier and 
Kostov in the analysis of the dense phase of polymers on random lattices 
of arbitrary genus [29]. They made use of the hermitian 0(n) multi-matrix 
model with a non-gaussian interaction, where 1/iV-expansion corresponds to 
the genus expansion and the n — > guarantees that a Hamiltonian cycle is 
connected [5]. In the present analysis of the compact polymer, I have written 
the generating function F{p,-^), whose 1/A^-expansion is again the genus 
expansion, in terms of the simplest matrix model: the hermitian gaussian 
1-matrix model without any limiting procedure. 
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5 Discussions 



I have obtained the generating function for the number of Hamiltonian cycles 
on surfaces of arbitrary genus. For genus g = 1, the contributions from 
winding and non- winding sector are determined. This can be done for an 
arbitrary genus by eq. (19). 

The hmit of large graphs (24), is consistent with the assertion that the 
system is in the same universahty class as the c = — 2 quantum gravity 
[5, 6]. In the way of calculating the generating function, I have made use of 
microscopic loop amplitudes of the gaussian hermitian 1-matrix model. 

Many deep connections have been found between the gaussian hermitian 
matrix model and the c = —2 quantum gravity since the first calculation of 
quantum gravity beyond the spherical limit at c = — 2 by Kostov and Mehta 
[31]. They have written the free energy at genus g in terms of correlation 
functions of the gaussian hermitian matrix model at order In contrast, 

F^^\p) in the present case gets contributions from correlators at all orders 
above The free energy in ref. [31] is the generating function for the 

number of maximal trees on graphs, while F^^^ (p) here generates the number 
of Hamiltonian cycles. 

It is interesting to compare eq.(24) with the number of connected trivalent 
fat graphs with n vertices in the absence of Hamiltonian cycles [5]. The latter 
behaves as [23, 32, 33]. 

^= E 1 X - (2 ■ 3^/^)'^nt^-i. (25) 

I introduce the 'random average' of ^{G) among G's in S"''^ by 

p(9) 

{n{G))GeSn,<, := 4^ ~ (2 • 3-3/4)n^(p+i)/2 ^ (0.877383 ••• )V^+'^/^ (26) 

This behavior may be compared to the case of the fixed fiat lattice with the 
same coordination number 3, i.e. the hexagonal lattice. The number Ti.{G) 
for the hexagonal lattice grows as [14, 15] 

n{G) - (33/V2)"n^ - (1.13975 • • • )"nT {n := #6"), (27) 

where the conformational exponent 7 is believed to be unity^. Having this 

^This definition of 7 is the standard one in the case of fixed lattice. One associate a 
base point to each Hamiltonian cycle in the present case, yielding an extra factor . 
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behavior of fixed regular lattice in mind, it is puzzling that the average 
{'H{G))GeS"'3 decreases as n grows. 

These totally different behaviors of (26) and (27) can be understood as 
follows. In the ensemble S*"'^ for arbitrarily large n, there are many G's 
with T-C{G) ~ 0. In fact, if a graph admits a Hamiltonian cycle, then the 
graph should be 2-edge-connected. Thus, as pointed out in ref. [5], G does 
not admit any Hamiltonian cycles if G is one-particle reducible. Therefore 
it is more natural to consider the random average over iS"'^, the restriction 
of 5*"'^ to one-particle irreducible graphs^. It is known that the number of 
one-particle irreducible graphs grows as [23] 

- E 1 X - (2-'/' • 3^/^)'^ni-i. (28) 

Therefore I obtain 

{'HiG))G^sn. (2'/' • 3-3/2)"n(^+^)/2 = (1.08866 • • • )"n(^+^)/l (29) 

It is interesting to note that the value 1.08866 ■ ■ ■ is near to the corre- 
sponding value 1.13975 • ■ ■ for the fixed hexagonal lattice as well as the field 
theoretic estimate 3/e = 1.10364- •• for regular lattices with coordination 
number three [34, 2]. For g — the power correction in (29) becomes n^, 
which coincides with that of the fixed flat hexagonal lattice (27). One is 
tempted to speculate that Hamiltonian cycles on random lattices resembles 
those on a regular lattice with the identical coordination number and the 
identical average curvature (zero for gf = 1) once one restricts oneself to 
one-particle irreducible graphs. 

In ref. [5], Eynard, Quitter and Kristjansen have mapped the problem of 
counting Hamiltonian cycles on planar random lattices to that of 0(n) model 
in the n — > limit on the sphere. Their Hamiltonian cycles are not furnished 
with directions and base points. The l/A/^-expansion of their model should 
contain the information of the number of Hamiltonian cycles on surfaces 
of arbitrary genus. I believe that their approach and the present one are 
complementary each other. 

Equation (22) formally resembles the expression for macroscopic loop 
amplitudes of the one-matrix model in ref. [35]. It may be interesting to 
translate the present analysis into the language of free fermions. 

^In this respect, I have aheady excluded disconnected graphs, which do not admit 
Hamiltonian cycles, in the definition of the random average (26). 
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